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Abstract 

Looking for an engine cycle with height output, multi-source of energy and less polluting pushes to reconsider the 
Stirling cycle. Several prototypes of engine were produced (Ford-Philips 4-215, Ross yoke, GPU-3... etc), but their 
performances remain weak compared with other types of internal combustion engine. In order to increase their 
performances and to analyze their operations, a numerical program of simulation taking into account thermal and 
mechanical losses was developed and a study of optimization of the design parameters was elaborated. The program 
which was applied to GPU-3 prototype of the General Motor gave results very close to the experimental results and 
leads to the optimization of the operating conditions. It also leads to the detennination of the optimal values of the 
geometrical and physical design parameters of the prototype and to the increase of its performances as long as the 
working gas pressure is maintained acceptable of the working gas in the engine. 
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1. Introduction 

The urgent need to preserve fossil fuels and to use renewable energies leads to using of the 
Stirling engines, these engines present an excellent theoretical output, equivalent to the Carnot one, they 
can use any source of energy (energy of combustion, solar energy...) and they are less polluting than the 
traditional engines (figure 1) (Stouffs P 1999) . 



Several prototypes were produced (Timoumi Y 2002), (Wang Jin T, 2002) and (S. Petrescu 2002) but 
their outputs remain very weak compared with the excellent theoretical yield. In fact, these engines have 
an extremely complex phenomena related to the compressible fluid mechanics, the thermodynamics and 
the heat transfer. The precise description and the comprehension of these highly non stationary 
phenomena are necessary in order to determine the different engine losses and to find the optimal values 
of the design variables. G. Popescu (1996) studied the influence of the heat losses and irreversibilities on 
the indices of the engine performance by a thermodynamic optimization which is endo- and exo- 
irreversible in finished time of the Stirling engine, showing that the most significant reduction in the 
performances is due to the no adiabatic regenerator. M. M. Salah El-Din (1999) and M. Costea (1999) 
studied the Stirling engines with solar energy; they concluded that the operating temperature and the 
engine performance depend on the internal irreversibility. I. Kolin (1991) showed that the discontinuous 
movement pushes the cycle to the constant volume lines. He studied several mechanisms of movement 
transformation and showed that more the movement is discontinuous more the real cycle approaches the 
ideal cycle. I. Tlili et al (2006) developed a dynamic model, taking into account the energy losses in the 
Stirling engine, the obtained results show that the losses by internal, external conduction and by pressure 
loss are most significant in the regenerator whereas the losses by Schuttle effect are more significant in 
the displacer piston. The other losses are negligible. The reduction of the losses increases the engine 
performance; it is necessarily based on the determination of the optimal geometrical and physical 
characteristics of the regenerator and the displacer piston. A study of optimization based on the dynamic 
model, I. Tlili et al (2006), is developed in this article. It allows to studying the influence of the 
geometrical and physical parameters on the prototype performance of Stirling engine and so to determine 
their optimal values. 
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2. Dynamic model of the Stirling engine 


The optimal geometrical and physical characteristics of the regenerator and the displacer piston. A study 
of optimization based on the dynamic model, I. Tlili et al (2006), is developed in this article. It allows to 
studying the influence of the geometrical and physical parameters on the prototype performance of 
Stirling engine and so to determine their optimal values. The dynamic model of the developed Stirling 
engine is based on the following assumptions: 

- The gas temperature in the different compartments of the engine is variable. 

- The temperatures of the cooler walls and the heater are taken constant. 

- The gas temperature in the different compartments is calculated referring to the law of perfect gases. 

The energy losses considered in this model are: 

- Energy Lost by internal conduction, through the solid matrix of the regenerator between the hot and cold 
parts which are calculated by the expression: 

A, 
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- Losses of energy by external conduction through the regenerator which is not adiabatic: 
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- The pressure losses by friction in the regenerator which leads to power losses: 
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- The losses of energy by Shuttle effect which comes from the shuttle movement of the displacer piston 
between the heater and the cooler: 
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Losses by irreversibility effect of the compression and the relaxation: 
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The dynamic model variables are given on the basis of the energy and mass conservation balance [11]: 
- The balance of energy conservation is written: 
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- The balance of the mass conservation in the engine gives: 

M=m ( j+m c +mf+m r +mj 1 (7) 

- By taking into account the losses by internal and external conduction in the regenerator, the quantities of 
heat exchanged in the regenerator are written: 

SWpcdr 
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- The work given by the cycle is: 
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3. Results of the model 

The developed model will be applied to a GPU-3 engine of the rhombic type, carried out by General 
Motor, Berchowitz D.M, Urieli 1 (1984). 

By introducing the geometrical and physical parameters into the models, the results are given in the 
following table 1: 

Table 1: Parameters and results of the prototype. 


prototype Parameters 

Values 

Regenerator thermal 
conductivity 

15 W/mK 

Regenerator porosity 

0.697 

Regenerator length 

0.0226 m 

Regenerator diameter 

0.0226 m 

Working fluid mass 

1.1362 g 

Displacer piston conductivity 

15 W/mK 

Displacer piston surface 

3.8 10' 3 m 2 

Displacer piston stroke 

0.046 m 

Calculated average pressure 

46.74 bar 

Calculated exchanged energy 

448.7206 J 

Calculated power 

4273 W 

Calculated efficiency 

39% 

Experimental power 

3500 W 

Experimental efficiency 

35 % 


The power and the efficiency calculated by the model are very close to the power and the real efficiency 
of the prototype. 

4. Optimization of the Stirling engine performances 

The energy losses are mainly observed in the regenerator, it represents 86 % of the total losses in the 
engine. They are primarily due to the losses by external and internal conduction and pressure losses. The 
energy lost by Schuttle effect in the displacer piston is also significant, it represents 13 %.The other losses 
are very weak [10]. The reduction of these losses generates an improvement of the engine 
performances .It is based on the determination of the optimal geometrical and physical parameters. 

These losses depend mainly on the matrix conductivity of the regenerator, its porosity, the in temperature 
variation, the working gas mass, the regenerator volume and the geometrical characteristics of the 
displacer. To study the influence of these parameters on the prototype performances, we change each time 
the studied parameter in the model and we maintain the other ones fixed, equal to the prototype 
parameters. 


4.1. Influence of parameters on the engine performances 

4.1.1. Effect of the conductivity> and the heat capacity’ of the regenerator matrix 


The performances of the engine vary according to the conductivity and to the heat capacity of material 
constituting the regenerator matrix. 
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Figure (2) shows that increase of matrix regenerator thermal conductivity leads to a reduction of the 
performances, which is explained by the increase of internal conduction losses in the regenerator 
according to conductivity I. Tlili et al (2006). 




Figure 2: Effect of regenerator thermal conductivity 
on the performances. 


Figure 3: Effect of the regenerator heat capacity on 
the perfonnances. 


Figure (3) shows that the engine performances are improved when the heat capacity of the regenerator 
matrix increases. 


The matrix of the regenerator can be made out of different materials. The performances of the engine are 
given according to the material nature in table 2. 


Table 2: Effect of the regenerator metal nature on the engine performances. 


Metal 

v cr 

(10 6 ) 

k 

W 

4 

Qr 

Steel 

3.8465 

46 

4258 

38.84 

441.25 

Stainless steel 

3.545 

15 

4273 

39.29 

448.72 

Copper 

3.3972 

389 

— 

— 

— 

Brass 

3.145 

100 

4080 

34.6 

415.67 

Aluminum 

2.322 

200 

3812 

29.16 

378.03 

Granite 

2.262 

2.5 

4091 

34.51 

430.75 

Concrete 

2.208 

0.92 

4081 

34.28 

429.9 

Glass 

2.1252 

1.2 

4062 

33.85 

427.8 

PVC 

1.45 

0.16 

3863 

30.04 

407.75 

Wood 

0.5125 

0.23 

3207 

23.75 

367.4 


Where: 

- V cr : Volumetric Capacity of heat J/m 3 K 

- k: Conductivity W/m.K 

- W : Engine power W 

- rp Engine efficiency % 

- Qr: Exchanged heat in the regenerator (j) 
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To increase heat exchange of the regenerator and to reduce the internal losses by conductivity, we should 
choose a material with high heat capacity and low conductivity. 

Steel, stainless steel, brass and granite are the best materials which can constitute the regenerator matrix. 

4.1.2. Effect of the fluid mass 

The increase of the fluid mass in the engine leads to a rise of the energy lost by pressure loss [13], 
however the engine output increases and the efficiency reaches a maximum which is about 40 % when the 
mass is equal to 0.8 g figure (4). 

When the mass increases, the efficiency reduction is due to the increase of the pressure loss and to the 
limitation of the heat exchange capacity in the regenerator and the heater. 

The use of a mass equal to 1.5 g in the engine allows having an acceptable output and a better power than 
in the prototype. 

4.1.3. Effect of the regenerator volume on the performances 

To vary the regenerator volume we can maintain the diameter fixed and vary the length or conversely. 
When the regenerator diameter is maintained constant equal to 0.0226 m, the length variation influences 
on the performances. The engine power and efficiency reach a maximum which corresponds to a length 
equal to 0.01 m then decrease quickly as shown in figure (5). That can be explained by the increase in the 
exchanged energy according to the length. 



Ms 

Figure 4: Effect of the fluid mass on the 
performances. 



Figure 5: Effect of the regenerator length on the 
engine perfonnances. 


When the regenerator length is constant L = 0.022 m. 


4.2. Optimization of the model parameters 


The Stirling engine performance depends on the choice of the geometrical and physical parameters. To 
study their influence on the engine performance, we gradually replaced these parameters by their optimal 
values. The obtained results are gathered in Table 3. 
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Table 3: Effect of the optimal values on the engine performances. 



optimal 

Value 

W 

9 

Qr 

Pmoy 

regenerator 

porosity 

65.5 % 

4554.2 

43.34 

474.31 

47.18 

regenerator 

length 

0.021 m 

4675.2 

41.04 

472.36 

47.93 

regenerator 

diameter 

0.024 m 

4765.1 

39.91 

475.76 

48.45 

working fluid 
mass 

1-15 g 

4849.6 

40.14 

480.71 

49 

displacer 

piston 

Conductivity 

1.2 

W/mK 

5079.4 

50.93 

479.87 

49.45 

displacer 
piston area 

3.86 

10' 3 m 2 

5183 

50.97 

483.71 

49.22 

displacer 
piston stroke 

0.047 m 

5106.2 

51.05 

472.67 

49.54 


Where: 

P moy : Average pressure (bar) 

At the beginning of table 3, we only replaced in the model, the prototype porosity (0.697) by the optimal 
porosity (0.655). The power and the efficiency are improved but the average pressure increases slightly. 
Therefore, we kept in the model the optimal porosity and we looked for the regenerator optimal length 
(0.021 m). By introducing it into the model, the power is improved and the efficiency remains acceptable. 
Then the optimal diameter (0.024 m) is used; the power and the exchanged energy in the regenerator 
increase but the efficiency decreases slightly. By replacing in the model the working fluid mass by the 
calculated optimal value (1.15g), the power and the efficiency increase but the average pressure remains 
acceptable. The displacer piston conductivity influences the loss by Shuttle effect, its reduction leads to a 
remarkable increase in the performances. 

By introducing in the model the optimal values of the area and the stroke of the displacer piston, the 
prototype performances are clearly improved. 

5. Conclusion 

The theoretical Stirling cycle presents an excellent theoretical output, but the realized prototypes have low 
outputs because of the considerable losses in the regenerator and the displacer piston, which are primarily 
due to the losses by external and internal conduction, the pressure losses in the regenerator and by 
Schuttle effect in the displacer piston. These losses depend on the geometrical and physical parameters of 
the prototype design. An optimization study of these parameters that was applied to GPU-3 engine of 
General Motor had a relatively significant real output (39%) and led to a reduction of these losses and to a 
notable improvement of the engine performance. To carry out this study, we introduced firstly the 
parameters of this prototype into a dynamic model; the obtained results were very close to the 
experimental results, Table 1. Then, we studied the influence of each geometrical and physical parameter 
on the engine performance and the exchanged energy of the regenerator. The reduction of the matrix 
porosity and conductivity of the regenerator increases the performance. The increase of the working fluid 
mass lead to an increase in the engine power and pressure but the efficiency reaches the maximum. When 
the displacer piston section increased and its stroke decreased, the engine output increases and the 
efficiency reaches the maximum. A low conductivity of the displacer piston reduces the losses by 
Schuttle effect and increases consequently the engine power and efficacity. Finally, we optimized these 
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parameters gradually by introducing them into the model and seeking at each time the optimum value. 
Although the real performance of the used prototype is relatively significant, the obtained results are 
considerably improved. The output passes from 39 % to 51 %, the power is improved of approximately 
20% and the average pressure is slightly increased as shown in Table 2. Applying this study to prototypes 
with weak performance or to new ones will lead to the determination of their design optimal parameters 
and consequently to interesting performance. 
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